THE STATISTICAL DISTRIBUTION OF THE ZEROS 
OF RANDOM PARAORTHOGONAL POLYNOMIALS 
ON THE UNIT CIRCLE 



MIHAI STOICIU 

Abstract. We consider polynomials on the unit circle defined by 
the recurrence relation 

$fc+i(z) = z$fe(z) - akKiz) fc > 0, $0-1 

For each n we take ao, ai, . . . , a„-2 i-i-d. random variables dis- 
tributed uniformly in a disk of radius r < 1 and q;„_i another 
random variable independent of the previous ones and distributed 
uniformly on the unit circle. The previous recurrence relation gives 
a sequence of random paraorthogonal polynomials {^n}n>o- Foi 
any n, the zeros of $„ are n random points on the unit circle. 

We prove that for any e*^ € the distribution of the zeros 
of in intervals of size O(^) near e'^ is the same as the distri- 
bution of n independent random points uniformly distributed on 
the unit circle (i.e., Poisson). This means that, for large n, there 
is no local correlation between the zeros of the considered random 
paraorthogonal polynomials. 



1. Introduction 

In this paper we study the statistical distribution of the zeros of 
paraorthogonal polynomials on the unit circle. In order to introduce 
and motivate these polynomials, we will first review a few aspects of 
the standard theory. Complete references for both the classical and the 
spectral theory of orthogonal polynomials on the unit circle are Simon 
[221 and [221. 

One of the central results in this theory is Verblunsky's theorem, 
which states that there is a one-one and onto map fi —>■ {an}n>o from 
the set of nontrivial (i.e., not supported on a finite set) probability 
measures on the unit circle and sequence of complex numbers {an}n>o 
with la^l < 1 for any n. The correspondence is given by the recurrence 
relation obeyed by orthogonal polynomials on the unit circle. Thus, 
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if we apply the Gram-Schmidt procedure to the sequence of polyno- 
mials l,z,z'^,... G L'^{dlD),dfi), the polynomials obtained ^o{z,dfj,), 
$1(2, dfi), $2 (-2, dfi) . . . obey the recurrence relation 

^k+iiz, dfi) = z^kiz, d/j,) - ak^liz, dfi) k>0 (1.1) 

where for ^k{z,dfi) = X]j=o^i'^"'' reversed polynomial $^ is given 
by ^l{z,dfi) = Yl^=o^k-jZ^ ■ The numbers from (jl.ip obey la^l < 1 
and, for any k, the zeros of the polynomial ^k+i{z,dfi) lie inside the 
unit disk. 

If, for a fixed n, we take ao, ai, ■ ■ ■ , an-2 ^ D and = P E d3 
and we use the recurrence relations (jl.l|) to define the polynomials 
$0) "^1) • • • ) ^n-i, then the zeros of the polynomial 

= z<^r.-iiz,dfi) -P^l_-^{z,dn) (1.2) 

are simple and situated on the unit circle. These polynomials (obtained 
by taking the last Verblunsky coefficient on the unit circle) are called 
paraorthogonal polynomials and were analyzed in jT2|, Jl]; see also 
Chapter 2 in Simon 

For any n, we will consider random Verblunsky coefficients by taking 
ao, «!, . . . , an-2 to be i.i.d. random variables distributed uniformly in 
a disk of fixed radius r < 1 and another random variable indepen- 
dent of the previous ones and distributed uniformly on the unit circle. 
Following the procedure mentioned before, we will get a sequence of 
random paraorthogonal polynomials {$„ = ^n{z,dfi, P)}n>o- For any 
n, the zeros of are n random points on the unit circle. Let us 
consider 

n 

k=l 

where zf^\ z^\ . . . , zi"^ are the zeros of the polynomial $„,. Let us also 
fix a point e*^ G D. We will prove that the distribution of the zeros 
of $„ on intervals of length 0{^) situated near e*^ is the same as the 
distribution of n independent random points uniformly distributed in 
the unit circle (i.e., Poisson). 

A collection of random points on the unit circle is sometimes called 
a point process on the unit circle. Therefore a reformulation of this 
problem can be: The limit of the sequence point process {Cn}n>o on a 
fine scale (of order 0(^)) near a point e*^ is a Poisson point process. 
This result is illustrated by the following generic plot of the zeros of 
random paraorthogonal polynomials: 
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This Mathematica plot represents the zeros of a paraorthogonal poly- 
nomial of degree 71 obtained by randomly taking ao, ai, . . . , aeg from 
the uniform distribution on the disk centered at the origin of radius | 
and random ajQ from the uniform distribution on the unit circle. On 
a fine scale we can observe some clumps, which suggests the Poisson 
distribution. 

Similar results appeared in the mathematical literature for the case 
of random Schrodinger operators; see Molchanov [201 Minami [T^. 
The study of the spectrum of random Schrodinger operators and of 
the distribution of the eigenvalues was initiated by the very important 
paper of Anderson |1] , who showed that certain random lattices exhibit 
absence of diffusion. Rigorous mathematical proofs of the Anderson 
localization were given by Goldsheid-Molchanov-Pastur ^1] for one- 
dimensional models and by Frohlich-Spencer JU] for multidimensional 
Schrodinger operators. Several other proofs, containing improvements 
and simplifications, were published later. We will only mention here 
Aizenman- Molchanov and Simon- Wolff |26^, which are relevant for 
our approach. In the case of the unit circle, similar localization results 
were obtained by Teplyaev ^7| and by Golinskii-Nevai [T^ . 

In addition to the phenomenon of localization, one can also ana- 
lyze the local structure of the spectrum. It turns out that there is 
no repulsion between the energy levels of the Schrodinger operator. 
This was shown by Molchanov [20] for a model of the one- dimensional 
Schrodinger operator studied by the Russian school. The case of the 
multidimensional discrete Schrodinger operator was analyzed by Mi- 
nami in |19j. In both cases the authors proved that the statistical 
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distribution of the eigenvalues converges locally to a stationary Pois- 
son point process. This means that there is no correlation between 
eigenvalues. 

We will prove a similar result on the unit circle. For any probability 
measure on the unit circle, we denote by {xo; Xi; X2, ■ ■ ■} the basis of 



dfi) obtained from {1, 



z, z 



^ , z 



. . .} by applying the Gram- 



Schmidt procedure. The matrix representation of the operator f{z) — > 
zf{z) on L^((9D, dfi) with respect to the basis {xo, Xi? X2, . . .} is a five- 
diagonal matrix of the form: 



/ "0 
Po 




V ... 



aipo 
— aiOo 

P2P1 




PiPo 
-Pitto 

-p2ai 







«3P2 P3P2 

-a3a2 -P3«2 



\ 



:i.4) 



(ao, ai, . . . are the Verblunsky coefficients associated to the measure /i, 
and for any n > 0, p„ = ^/l^-\oin^). This matrix representation is a 
recent discovery of Cantero, Moral, and Velasquez |H] . The matrix C is 
called the CMV matrix and will be used in the study of the distribution 
of the zeros of the paraorthogonal polynomials. 

Notice that if one of the a's is of absolute value 1, then the Gram- 
Schmidt process ends and the CMV matrix decouples. In our case, 
|q;„_i| = 1, so pn-i = and therefore the CMV matrix decouples be- 
tween {n — 1) and n and the upper left corner is an (n x n) unitary 
matrix C*-"\ The advantage of considering this matrix is that the zeros 
of are exactly the eigenvalues of the matrix C*-"^ (see, e.g., Simon 
[221). ^ill some techniques from the spectral theory of the 
discrete Schrodinger operators to study the distribution of these eigen- 
values, especially ideas and methods developed in [T], j2], 0, [HI, f^ . 
[21]. However, our model on the unit circle has many different features 
compared to the discrete Schrodinger operator (perhaps the most im- 
portant one is that we have to consider unitary operators on the unit 
circle instead of self-adjoint operators on the real line). Therefore, we 
will have to use new ideas and techniques that work for this situation. 

The final goal is the following: 

Theorem 1.1. Consider the random polynomials on the unit circle 
given by the following recurrence relations: 



> 0, 



$0 = 1 



(1.5) 
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where ao, ai, ■ ■ ■ , an~2 are i.i.d. random variables distributed uniformly 
in a disk of radius r < 1 and a;„__i is another random variable indepen- 
dent of the previous ones and uniformly distributed on the unit circle. 

Consider the space Q = {a = (oo, ai, . . . , an-2, «n-i) ^ -D(0,r) x 
D{0,r) X ■ ■ ■ X D{0,r) X dll>} with the probability measure P obtained by 
taking the product of the uniform (Lebesgue) measures on each D{0,r) 
and on dB>. Fix a point e''^° G D and let (^("^ be the point process defined 
by (113). 

Then, on a fine scale [of order ^) near e^^^ , the point process (^"^^ con- 
verges to the Poisson point process with intensity measure n |^ {where 
1^ is the normalized Lebesgue measure). This means that for any fixed 
ai < hi < a2 < b2 < ■ ■ ■ < am < bjji and any nonnegative integers 
ki,k2, . . . , km, we have: 

(1.6) 

as oo. 

2. Outline of the Proof 

From now on we will work under the hypotheses of Theorem 11.11 We 
will study the statistical distribution of the eigenvalues of the random 
CMV matrices 

C(^)=C^^) (2.1) 

for a G (with the space Q defined in Theorem II. Ij) . 

A first step in the study of the spectrum of random CMV matrix is 
proving the exponential decay of the fractional moments of the resol- 
vent of the CMV matrix. These ideas were developed in the case of 
Anderson models by Aizenman-Molchanov |T| and by Aizenman et al. 
[2] . In the case of Anderson models, they provide a powerful method for 
proving spectral localization, dynamical localization, and the absence 
of level repulsion. 

Before we state the Aizenman-Molchanov bounds, we have to make 
a few remarks on the boundary behavior of the matrix elements of the 
resolvent of the CMV matrix. For any 2; G D and any < k,l < (n — l), 
we will use the following notation: 

(2.2) 

kl 



FkAz.c^:^] 



d"^ + z 



in) 
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As we will see in the next section, using properties of Caratheodory 
functions, we will get that for any a G fi, the radial limit 

Fki{e'\C^^^) = limFH(re*^C(")) (2.3) 
exists for Lebesgue almost every 

for any s G (0, 1). Since the distributions of ao, ai, . . . , On-i are ro- 
tationally invariant, we obtain that for any fixed e*^ G 9©, the radial 
limit Fki{e^^ ,C^^) exists for almost every a G fi. We can also define 



(2.4) 

kl 



and 



GfcKe'^Ci")) = limGfcKre^^Ci")) (2.5) 

rjl 



Using the previous notation we have: 

Theorem 2.1 (Aizenman-Molchanov Bounds for the Resolvent of the 
CMV Matrix) . For the model considered in Theorem li.il and for any 
s G (0, 1), there exist constants Ci,Di > such that for any n > 0, 
any k,l, < k,l < n — 1 and any e*^ G 83, we have: 

E(|FfeKe'',Ci"))f) <Cie-^il^-'l (2.6) 

where C*^"^ is the {n x n) CMV matrix obtained for ao, ai, . . . q;„_2 uni- 
formly distributed in D{0,r) and uniformly distributed in 83. 

Using Theorem 12.11 we will then be able to control the structure of 
the eigenfunctions of the matrix C^^\ 

Theorem 2.2 (The Localized Structure of the Eigenfunctions). For 

the model considered in Theorem M.R the eigenfunctions of the random 
matrices C-"^ = Ci"^ are exponentially localized with probability 1, that 
is exponentially small outside sets of size proportional to (ln?7,). This 
means that there exists a constant D2 > and for almost every a G fi, 
there exists a constant Ca > such that for any unitary eigenfunction 
tf^a \ there exists a point m{(p^'') (1 < m{if^a '') < n) with the property 
that for any m, \m — m{(p^a'')\ > D2 ln{n + 1), we have 

|(/?i")(m)| < C7^e-('/^^)l'"-'"(^^"')l (2.7) 
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The point m{ip')^^) will be taken to be the smallest integer where the 
eigenfunction (^^"^(m) attains its maximum. 

In order to obtain a Poisson distribution in the limit as n ^ oo, we 
will use the approach of Molchanov |j20j and Minami jTH]- The first 
step is to decouple the point process C^""^ into the direct sum of smaller 
point processes. We will do the decoupling process in the following way: 
For any positive integer ra, let C-"-* be the CMV matrix obtained for 
the coefficients Q^i, • • • , with the additional restrictions ar^i = 

[ In n J 

e*''!, a^r n 1 = e*''^, = e*''''""i where e*''\ e*''^, . . . , e*''['""i are 

independent random points uniformly distributed on the unit circle. 
Note that the matrix C*^"-* decouples into the direct sum of ~ [Inn] 
unitary matrices Cj"'', Cg""*, • • • , C^^^l^y We should note here that the 

actual number of blocks C^'' is slightly larger than [In n] and that the 
dimension of one of the blocks (the last one) could be smaller than 

[inn] • 

However, since we are only interested in the asymptotic behavior 
of the distribution of the eigenvalues, we can, without loss of gener- 
ality, work with matrices of size = [Inn] [j^] . The matrix C'-^-* is 

the direct sum of exactly [Inn] smaller blocks c'i^\ C^i^l^y 

We denote by C^^'^^ = E£=^"^ '^.(p) where z^^\ z!f\..., z[^]^^^^ are the 

eigenvalues of the matrix Cp^\ The decoupling result is formulated in 
the following theorem: 

Theorem 2.3 (Decoupling the point process). The point process C*-^-* 
can be asymptotically approximated by the direct sum of point processes 
Yjp=i C^^'^'^ ■ other words, the distribution of the eigenvalues of the 
matrix C^^^ can be asymptotically approximated by the distribution of 
the eigenvalues of the direct sum of the matrices C[^^ , C'^'^ , ■ ■ ■ , C^^^j . 

The decoupling property is the first step in proving that the statis- 
tical distribution of the eigenvalues of C^^^ is Poisson. In the theory 
of point processes (see, e.g., Daley and Vere- Jones |7|), a point process 
obeying this decoupling property is called an infinitely divisible point 
process. In order to show that this distribution is Poisson on a scale of 
order O(-) near a point e*^, we need to check two conditions: 



[Inn] 



P(^Wp)(^(]V,0)) > 1) ^ |A| as n^oo (2. 
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[Inn] 

ii) ^P(C(^'^H^(^,^)) >2) ^0 as n^oo (2.9) 
p=i 

where for an interval A = [a, b] we denote by A{N, 6) = (e*(^+^\ e*'-^^^^) 
and I ■ I is the Lebesgue measure (and we extend this definition to unions 
of intervals). The second condition shows that it is asymptotically im- 
possible that any of the matrices Cj^-* , C'^^ , • • • , C^^^j has two or more 
eigenvalues situated an interval of size Therefore, each of the ma- 
trices C'^\ C'^\ . . . , Cp^^j contributes with at most one eigenvalue in 

an interval of size But the matrices C'^\ • • • , C^^^j are decou- 
pled, hence independent, and therefore we get a Poisson distribution. 
The condition i) now gives Theorem 11.11 

The next four sections will contain the detailed proofs of these the- 
orems. 



3. AlZENMAN-MOLCHANOV BOUNDS FOR THE RESOLVENT OF THE 

CMV Matrix 

We will study the random CMV matrices defined in (j^.lj) . We will 
analyze the matrix elements of the resolvent (C*^"-* ~ ^)~^ of CMV 
matrix, or, what is equivalent, the matrix elements of 

F{z, C(")) = (C(") + z)(C(") - z)-^ = I + 2z (C(") - z)-^ (3.1) 

(we consider z G D). More precisely, we will be interested in the expec- 
tations of the fractional moments of matrix elements of the resolvent. 
This method (sometimes called the fractional moments method) is use- 
ful in the study of the eigenvalues and of the eigenf unctions and was 
introduced by Aizenman and Molchanov in jT]. 

We will prove that the expected value of the fractional moment of 
the matrix elements of the resolvent decays exponentially (see (j2.6p ). 
The proof of this result is rather involved; the main steps will be: 

Step 1. The fractional moments E 

bounded f Lemma l3.1|) . 



Step 2. The fractional moments E 
uniformly along the rows f Lemma 13. 6p . 

Step 3. The fractional moments E 
tially (Theorem 1211). 

We will now begin the analysis of E 



kl 



Zj Ca 



are uniformly 
converge to 
decay exponen- 
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It is not hard to see that Re [(C^"^ + z)(C^"-^ — z)~^~\ is a positive 
operator. This will help us prove: 

Lemma 3.1. For any s G (0,1), any k,l, 1 < k,l < n, and any 

z e D U we have 

E(|F,K^,C("))f) <C (3.2) 

where C = L . 

Proof. Let F^{z) = {ip, (Ci"^ + z)(d"^ - z)-^(p). Since ReF^ > 0, the 
function is a Caratheodory function for any unit vector ip. Fix 
p G (0, 1). Then, by a version of Kolmogorov's theorem (see Duren jH] 
or Khodakovsky T7]), 

'^|(^,(Ci")+pe^')(Ci")-pe^')"V)rf < C, (3.3) 



where Ci = 

The polarization identity gives (assuming that our scalar product is 
antilinear in the first variable and linear in the second variable) 

1 ^ 

F,Kpe^^C(")) = - $^(-0™ ((5. + ^"^50, F(pe^^C("))(5, + ^50) 

m=0 

(3.4) 

which, using the fact that |a + 6|'*<|a|* + |6|'', implies 



m=0 \ V V 

(3.5) 

Using (|3.3|) and (|3.5|) . we get, for any ci"\ 

^'>.Kpe^C("))r|<C (3.6) 



where C = -^-^s- 

Therefore, after taking expectations and using Fubini's theorem, 

^"E(|FH(pe^^CW)r)^ <C (3.7) 

The coefficients ao, ai, . . . , define a measure dfi on 5D. Let us 
consider another measure (i/ie(e*'^) = (i/i(e**^'^~^^). This measure defines 
Verblunsky coefficients ao.e, «i,6», • • • , ctn-i.e, a CMV matrix C^"^, and 
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unitary orthogonal polynomials ipQ^, ipi^e . . . , fn-i,9- Using the results 
presented in Simon ^22^, for any k, < k < n — 1, 



ak,6 



(3.8) 
(3.9) 



The relation ()3.9j) shows that for any k and 9, Xk,e{z) = Xkfi Xfc(e "^^z) 
where \\kfi\ = 1. 

Since ao, ai, . . . , are independent and the distribution of each 
one of them is rotationally invariant, we have 



E(|Ffe,(pe^C("))f) =E 
But, using dnUI) and ^U^ . 



Fkiipe'^XS' 



(3.10) 



Fki{pe'',Ci: 



_ p^,e Ueien XkAe^n dpoien 



XiAenXkAe'^dKe'^^-'^) 



+ pe 



ie 



dl 



(,i{T+e) _ pg 

e*" + p 



Xi,e{e'^'^'^)Xk,e{e'^^^'^)dp{e 



e" — p 



Xiien Xkie^^dpie^n 



where |A;,eAfc,e| = 1- 

Therefore the function 6 
ing ()3.7|) . we get 

E 



E 



FM(pe^^a' 



IS constant, so, us- 



< C 



(3.11) 



Since p and 6 are arbitrary, we now get the desired conclusion for any 

zeB). 

Observe that, by ()3.4j) . Fki is a linear combination of Caratheodory 
functions. By ^8j, any Caratheodory function is in if'^(©) (0 < s < 1) 
and therefore it has boundary values almost everywhere on 83. Thus 
we get that, for any fixed a E Q and for Lebesgue almost any z = e^^ E 
dJ}, the radial limit Fki{e^^ ,C^^) exists, where 



F,Ke^^C^"^) = limFH(pe^^C(")) 
pTi 



(3.12) 
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Also, by the properties of Hardy spaces, ■ , Ci"'*) G L^{d3) for 
any s G (0, 1). Since the distributions of ao, ai, . . . , On-i are rotation- 
ally invariant, we obtain that for any fixed e*^ G 9D, the radial limit 
-Ffc/(e*^, Ci"''*) exists for almost every a G fi. 

The relation ()3.11|) gives 




(3.13) 



By taking p ] 1 and using Fatou's lemma we get: 

E(|F,Ke^',Ci"))f) <C (3.14) 

□ 

Note that the argument from Lemma lTTl works in the same way when 
we replace the unitary matrix ci"^ with the unitary operator Ca (cor- 
responding to random Verblunsky coefficients uniformly distributed in 
D{0,r)), so we also have 

E(|FH(e*^On <C (3.15) 

for any nonnegative integers k, I and for any e*^ G dB>. 

The next step is to prove that the expectations of the fractional 
moments of the resolvent of C*^"-* tend to zero on the rows. We will 
start with the following lemma suggested to us by Aizenman jS]: 

Lemma 3.2. Let {Xn = Xn{uj)}n>o, uj E Q be a family of positive 
random variables such that there exists a constant C > such that 
E(X„) < C and, for almost any uj E Q, lim„,_^oo -^n(i^) = 0. Then, for 
any s G (0, 1), 

lim E(X^) = (3.16) 

n— >oo 

Proof Let e > and let M > such that M^~^ < e. Observe that if 
Xn{uj) > M, then X^(cu) < M'-^Xn{uj). Therefore 

X^^ico) < X^^ico) X{.;X„H<Af}(c^) + M^-iX„(cu) (3.17) 

Clearly, E(M*^^X„) < eC and, using dominated convergence, 

E(X^XV;X„H<M}) ^ as n-^oo (3.18) 

We immediately get that for any £ > we have 

limsupE(X^) < E(X:x{.;X„H<A/}) +£C (3.19) 

n— >oo 

SO we can conclude that ()3.16p holds. □ 
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We will use Lemma ltT^ to prove that for any fixed J, E ( | Ca) \ '*) 

and E (|Fjj+fc(e*^, Cq) I'') converge to as A; ^ cxd. From now on, it 
will be more convenient to work with the resolvent G instead of the 
Caratheodory function F. 

Lemma 3.3. Let C = be the random CMV matrix associated to a 
family of Verblunsky coefficients {an}n>o with an i.i.d. random vari- 
ables uniformly distributed in a disk D{0,r), < r < 1. Let s G (0, 1), 
2; e D U 83, and j a positive integer. Then we have 

limE(|G',-,+fc(^,C)r) = (3.20) 

Proof. For any fixed ^ G D, the rows and columns of G{z,C) are 
at infinity, hence converge to 0. Let s' G (s,l). Then we get ()3.20p 
applying Lemma [3.21 to the random variables Xk = \Gjj+k{z,C)f and 
using the power 7 < 1. 

We will now prove ()3.20j) for z = e*^ G In order to do this, 

we will have to apply the heavy machinery of transfer matrices and 
Lyapunov exponents developed in [23]. Thus, the transfer matrices 
corresponding to the CMV matrix are 

Tr,{z) = A{an,z)...A{ao,z) (3.21) 

where A{a, z) = {1 — lap)^"*^/^ ( ~j^") and the Lyapunov exponent is 

7(^)= hm -log||T„(z,K})|| (3.22) 

n— >oo fl 

(provided this limit exists). 

Observe that the common distribution d^a of the Verblunsky coeffi- 
cients a„ is rotationally invariant and 

/ — log(l — u;) (i/iQ,(u;) < 00 (3.23) 
Jd(o.i) 



and 



/ — log Ic^l (i/ia(i^) < C)0 (3.24) 

Jd{0,1) 

Let us denote by du^ the density of eigenvalues measure and let Lf^"^ 
be the logarithmic potential of the measure du^, defined by 

^'^"(e'')=/ log I ^ , . du^ien (3.25) 
J do \e'^-e"\ 

By rotation invariance, we have dvjq = ^ and therefore U'^'^'^ is 
identically zero. Using results from ,2^ , the Lyapunov exponent exists 
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for every z = e*^ G dV> and the Thouless formula gives 

7(^) = -i / \og{l-\uj\^)dfi^{u) (3.26) 
Jd(o,i) 



By an immediate computation we get 7(2) = 2r°^^^~^ ^ ^■ 

The positivity of the Lyapunov exponent 7(e*^) imphes (using the 
Ruelle-Osceledec theorem; see [23^) that there exists a constant A 7^ 1 
(defining a boundary condition) for which 



hmT„(e*^)( =0 (3.27) 

From here we immediately get (using the theory of subordinate so- 
lutions developed in |2S1) that for any j and almost every e*^ G d3, 

lim G'jj+fc(e*^C) = (3.28) 

We can use now ()3.15|) and ()3.28p to verify the hypothesis of Lemma 
for the random variables 



Xk=\G,,,+k{e'',C)\' (3.29) 
where s' G (s, 1). We therefore get 

lim E (IG. ,+fc(e^^ C)\") = (3.30) 

fc— >oo 

□ 

The next step is to get the same result for the finite volume case 
(i.e., when we replace the matrix C = Cq by the matrix C^^). 

Lemma 3.4. For any fixed j , any s G (0, |), and any z G © U dlD), 

lim E(|G,,+.(^,C("))n = (3.31) 

Proof. Let C be the CMV matrix corresponding to a family of Verblun- 
sky coefficients {an\n>o, with |a„| < r for any n. Since E {\Gjjj^k{,z, C)|^) 
and E (|Gjj+fc(2;,C)|^'') as A; — 00, we can take fc^ > such that 
for any A; > A;^, E {\Gj^j+k{z,C)\') < e and E (|Gjj+fc(z, C)|^") < e. 

For n > (fee + 2), let C*^"^ be the CMV matrix obtained with the same 
ao, a„_2, a„, . . . and with a„_i G d3. From now on we will use 

G{z,C) = {C- z)-^ and G{zXa^) = {Ct^ - z)~\ Then 

(C(") - ;2)'i - (C - z)-i = (C - z)-i(C - Ci"))(Ci") - ^)-i (3.32) 
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Note that the matrix (C — C''"^) has at most eight nonzero terms, 
each of absolute value at most 2. These nonzero terms are situated at 
positions (m, m') and I'm — n\ < 2, |m' — n\ < 2. Then 

8 terms 

(3.33) 

Using Schwartz's inequality, 

E(ic-.)-^n(e)-.)„}^.^,r)< 



We clearly have m > k,, and therefore E (|C — -z)^^^'^) — ^- Also, 
from Lemma l3.H there exists a constant C depending only on s such 

that E (id") - < C. 

Therefore, for any > fc„ E (|(d") - < £ + e^^^C. 

Since e is arbitrary, we obtain ()3.3H) . □ 

Note that Lemma 13.41 holds for any s G (0, |). The result can be 
improved using a standard method: 

Lemma 3.5. For any fixed j , any s G (0, 1), and any z eB), 

hm E(|G„+.(z,C("))|^) = (3.35) 

k^oo, k<n 

Proof. Let s G [i,l),t G (s, l),r G (0, |). Then using the Holder 
inequality for p = and for q = j^, we get 



(3.36) 



From Lemma EH E( | (Ci^^ — ^)jj^+fc|*) is bounded by a constant de- 
pending only on t and from Lemma IS31 E(|(Ci"') — z)jj^j^\^) tends to 
as ^ oo. We immediately get (|3.35|) . □ 

We can improve the previous lemma to get that the convergence to 
of E(|(d") 

— ^)jj^_i_fe|^) is uniform in row j. 
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Lemma 3.6. For any e > 0, there exists a > such that, for 
any s, k,j, n, s E (0, 1), k > k^,n > 0,0 < j < {n — 1), and for any 
^ e D U (9D, we have 



E G 



< e 



(3.37) 



Proof. As in the previous lemma, it is enough to prove the result for 
all 2; G D. Suppose the matrix C*^"^ is obtained from the Verblunsky 
coefficients ao,ai, . . . ,an-i- Let's consider the matrix C^l obtained 
from the same Verblunsky coefficients with the additional restriction 
ttm = e*^ where m is chosen to be bigger but close to j (for example 
m = j + 3). We will now compare (C^") - -2)"]+^ and {C^l - -2)^)+^- 
By the resolvent identity, 



+k\ 



< 2 



jj+k 

E 

\l-m\<2,\l'-m\<2 



J,3' 



(3.38) 



(3.39) 



The matrix (C^"^ — z)^^ decouples between m — 1 and m. Also, since 
\V — m\ < 2, we get that for any fixed e > 0, we can pick a ke such that 
for any k > k^ and any \l' — m| < 2, we have 



< e 



(3.40) 



(In other words, the decay is uniform on the 5 rows m — 2,m — l,m,m + 
1, and m + 2 situated at distance at most 2 from the place where the 
matrix C^l decouples.) 

As in Lemma (3.41 we can now use Schwartz's inequality to get that 
for any e > and for any s G (0, |) there exists a k^ such that for any 
j and any k > k;,. 



e(|(c 



'j,'j+k\ 



< £ 



(3.41) 



Using the same method as in Lemma 13.51 we get ()3.37|) for any 
SG (0,1). 

□ 



We are heading towards proving the exponential decay of the frac- 
tional moments of the matrix elements of the resolvent of the CMV 
matrix. We will first prove a lemma about the behavior of the entries 
in the resolvent of the CMV matrix. 
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Lemma 3.7. Suppose the random CMV matrix C^"^^ = is given as 
before {i.e., ao, «i, . . . , an_2, are independent random variables, 
the first {n — 1) uniformly distributed inside a disk of radius r and the 
last one uniformly distributed on the unit circle) . Then, for any point 
e*^ G and for any a G where G{e^^ ,C^^) = {C^^ — e^^)^^ exists, 
we have 



n (pie r^^-h 



\k-i\ + \l-i\ 

< ( ) (3-42) 



r 



Proof. Using the results from Chapter 4 in Simon |22], the matrix el- 
ements of the resolvent of the CMV matrix are given by the following 
formulae: 

\(r-y^-^^ = ! ('^^y^Xiiz)pkiz), k>l or k = l = 2n-l 
^ ^ki \ {2z)-^7ii{z)xkiz), l>k or k = l = 2n 

(3.43) 

where the polynomials Xi{z) are obtained by the Gram-Schmidt pro- 
cess applied to {1, z, z~^, . . .} in L^(c}D, dfi) and the polynomials Xk{z) 
are obtained by the Gram-Schmidt process applied to {1, z~^ , z . . .} in 
L'^{dlD),dfj,). Also, pn and 7r„ are the analogs of the Weyl solutions of 
Golinskii-Nevai ^3] and are defined by 

Pn = yn + F{z)Xn (3.44) 

7r„ = T„ + F{z)xn (3.45) 

where ?/„ and T„ are the second kind analogs of the CMV bases and 
are given by 

/ z"V2i n = 2l , . 

= [ -z-^r^i-, n = 2l-l ^3.46) 

\ z-'+'^2i-i n = 2l-l ^'^■^'^ 

The functions ipn are the second kind polynomials associated to the 
measure /i and F{z) is the Caratheodory function corresponding to /i 
(see PI). 

We will be interested in the values of the resolvent on the unit circle 
(we know they exist a.e. for the random matrices considered here). 
For any z G 8^, the values of F{z) are purely imaginary and also 

Xn{z) = Xn{z) and Tn{z) = -yn{z). In particular, \Xn{z)\ = \xn{z)\ 
for any z G 83. 

Therefore Uniz) = Tn{z) + F{z)xniz) = -Pn{z), so \-Kn{z)\ = \Pn{z)\ 
for any z G dB. We will also use |x2n+i(^)| = \'f2n+iiz)\,\x2niz)\ = 
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\'f*2niz)\Ax2n{z)\ = \'f2niz)\, and \x2n-iiz)\ = Iv^L-il^^)! ^oi any z e 
d3. Also, from Section 1.5 in [22], we have 



Vn[,Z) 



< c 



(3.48) 



for any z G 9D, where C = 2/vl 

The key fact for proving ()3.48jl is that the orthogonal polynomials 
ipn satisfy a recurrence relation 

fn+l{z) = p-\zifn{z) - anfliz)) (3.49) 

This immediately gives the corresponding recurrence relation for the 
second order polynomials 

Ipn+liz) = p~^{zilJn{z) + an-^niz)) (3.50) 

Using (j3.49|) and (j3.5(Jj) . we will now prove a similar recurrence rela- 
tion for the polynomials vr^. For any z G dH), we have 

■^21+iiz) = T2i+iiz) + F{z)x2i+iiz) 

= Z-\ij2l+liz) + F{z)if2l+l{z)) 



-p2i^ Zn2l{z) +P2i^ a2lTT2l{z) 



and similarly we get 



Tr2l{z) = -p2i_i7l2l-l{z) - a2l-lp2l-iTr2l-l{z) 



(3.51) 



(3.52) 



where we used the fact that for any 2; G D, F{z) is purely imaginary, 
hence F{z) = —F{z). 

Since Pn^< :^7f=f ^ the equations ()3.5H) and ()3.52p will give that for 
any integer n and any z G ©, 

T^n±l{z) 



vr„ U 



< c 



(3.53) 



where C = 2/ \/l — r'^. 

Using these observations and ()3.43|) we get, for any z G (9D, 



[(c(") - zr] 



< c 



[(cw - zr] 



k,l±l 



and also 



We can now combine ()3.54|) and ()3.55j) to get ()3.42p . 



(3.54) 

(3.55) 

□ 



We will now prove a simple lemma which will be useful in computa- 
tions. 
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Lemma 3.8. For any s G (0, 1) and any constant P E C, we have 

1 1 



Proof. Let (3 = (3i + i(32 with /32 G M. Then 



dx < — — dx (3.56) 



dx = / — tt; rrr^ dx < — dx 



-1 k - J -I - Pi)' + Pl\'^' ~ J -I k - A 

(3.57) 

But l/|a;|* is the symmetric decreasing rearrangement of l/\x — (3iY 
so we get 

1 I 

dx< T^dx (3.58) 



|x /^i I J —1 1*^1 
and therefore we immediately obtain ()3.56p . □ 

The following lemma shows that we can control conditional expec- 
tations of the diagonal elements of the matrix C^"^. 

Lemma 3.9. For any s G (0, 1), any k, 1 < k < n, and any choice of 
ao, ai, . . . , ajt-i, ak+i, . . . , a;„_2, a;„_i, 

E(|F,fc(z,C("))f I <C (3.59) 

where a possible value for the constant is C = 32*. 

Proof. For a fixed family of Verblunsky coefficients the diag- 

onal elements of the resolvent of the CMV matrix C can be obtained 
using the formula: 

{6k,{C + z){C-z)-'6,)= [ 4^l^fc(e^')r^^/i(e^') (3.60) 



where /i is the measure on dV> associated with the Verblunsky coeffi- 
cients {anjneN and {v?n}neN are the corresponding normalized orthog- 
onal polynomials. 



Using the results of Khrushchev |THj , the Schur function of the mea- 
sure |v9fc(e*^)p (i/x(e*^) is: 

gk{z) = f{z; ak, ak+i, . . .) f{z; -Ok-i, -ak-2, -^o, 1) (3.61) 

where by f{z; S) we denote the Schur function associated to the family 
of Verblunsky coefficients S. 

Since the dependence of f{z; ak, dk+i, • • •) on ak is given by 

N ak + zf{z;ak+i,ak+2---) 

f{z;ak,ak+i,...) = —^ — r (3.62) 

1 + akzf{z; ak+i,ak+2 ■ ■ ■) 
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we get that the dependence of gk{z) on ak is given by 

9k{z) 

where 

Cl = f{z] -ttfc-l, -ayk-2, • • • , -tto, 1) 



^ 1 + 



(3.63) 



C2 = zf{z; afc+i,afc+2, 



(3.64) 
(3.65) 



Note that the numbers Ci and C2 do not depend on a^, \Ci\, IC2I < 1. 

We now evaluate the Caratheodory function F{z; |y9jt(e*^)p (i/i(e*^)) 
associated to the measure |v?/c(e*^)p (i/i(e*^). By definition, 

e'^ + z 



Fiz;\Men\'dfiie^')) 

We now have 
\F{z-\Me'')\'dfi{e^'))\ = 

It suffices to prove 

sup 

u 

Clearly 



{5k,iC + z)iC-z)-'5,) 



1 + zgk{z) 



1 - zgk{z) 



< 



[)i,«J2eD JD(0,r) 



1 _ ^. ^ 



(3.66) 
(3.67) 

(3.68) 
(3.69) 




. - / X (3.70) 

1 + akW2 - wi{ak + W2) 

For ak = x + iy, 1 + 0^^2 - Wi{ak + W2) = x{-Wi + W2) +y{-iwi - 
iwi) + (1 — W1W2). Since for wi,W2 € D, {—Wi + W2), {—iwi — 1102), 
and (1 — W1W2) cannot be all small, we will be able to prove ()3.69|1 . 

If I - wi + W2I > e, 



I 

J D{Q,r) 



l+a 



kW2 



dai. < 



1 



\x + yD + E\ 



■ dx dy 



(3.71) 



< 2 



dx 



\x\ 



1 — s \e 



(3.72) 
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(where for the last inequality we used Lemma f3.8|l . 
The same bound can be obtained for \wi + W2I > s. 
If I — Wi + W2I £ and \wi + W2I ^ then 

\x{-wi + W2) + y{-iwi - iw2) + (1 - ^1^2) I > (1 - £^ - 4e) (3.73) 
so 



dak < 2'+^ 



1 



Therefore for any small e, we get (|3.59|) with 



C = max 



1 — s V £ 



For example, for e = 1/8, we get C = 32'^. 
We will now be able to prove Theorem 12.11 



(3.74) 



(3.75) 



□ 



Proof of Theorem \2. 11 We will use the method developed by Aizenman 
et al. for Schrodinger operators. The basic idea is to use the uni- 
form decay of the expectations of the fractional moments of the matrix 
elements of C^"'^ (Lemma I3.(i|l to derive the exponential decay. 

We consider the matrix C'-"^ obtained for the Verblunsky coefficients 
ao, ai, . . . , Fix a k, with < k < {n — 1). Let C["'^ be the 

matrix obtained for the Verblunsky coefficients a;o,«i, . . . ,«n-i with 
the additional condition ctfc+m = 1 and C^-' the matrix obtained from 
ao, an-i with the additional restriction ak+m+s = e*^ ("^ is an 
integer > 3 which will be specified later, and e*^ is a random point 
uniformly distributed on 83). 

Using the resolvent identity, we have 

(CW-z)-i-(ci")-z)-i = (cJ")-z)-i(c;"^-C("))(C(")-z)-i (3.76) 
and 

(CW-^)-i-(C^")-^)-i = (C(")-^)-i(Cf)-C("))(C^")-z)-i (3.77) 
Combining (j3.76j) and (j3.77|) . we get 

(C(") - zy' = (cf") - z)-' + (cj") - z)-' (ci"^ - c(")) (c^") - z)-' 

+ (Ci") - z)-' (Ci") - C(")) (C(") - z)~^ (C^"^ - C(")) (4") - z)-' 

(3.78) 
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For any k, I with / > (A; + m), we have 



in) ^ .-I 







J kl 



and 



kl 



(3.79) 
(3.80) 



(C;"^-^)-^(Cj")-C("))(C?^-^)-i 

Therefore, since each of the matrices {C[^^ — C^"-*) and (C2"'' — C) has 
at most eight nonzero entries, we get that 

[(c'"'-^)-'].,= 

E <' - Cl"' - (c'"' - .-);,■,. (4"' - c<"'),„. (c<»> - 

64 terms 

(3.81) 

which gives 

e(|(cW-^),/ 



<4^ 5^ E(|(Cf) 



Iks 



1 //7(ri.) 



-1 //^(n) 



(3.82) 



64 terms 

where, since the matrix Cj"^ decouples at {k + m), we have |s2 — (fc + 
m)| < 2 and, since the matrix C["^ decouples at (A; + m + 3), we have 
|S3 - (A; + m + 3)| < 2. 

By Lemma (3. 7( we have for any e*^ G 83, 



(C(") 



< 



/fc+m+l,fc+m+l I 



(3.83) 



Observe that {€[ — z)^/^ and (C2 — z)^^"^ do not depend on ak+m+i, 
and therefore using Lemma f3. 91 we get 



E 



< 



32' 



(n) 



Since the random variables (Cj"'* — z)^,^'^ and (Cg 



,{n) 



(3.84) 



z)^^] are inde- 



pendent (they depend on different sets of Verblunsky coefficients), we 
get 



E 



S41 



< C(.,r)E (|(Cj") - .),-\[) E ([(Ci") - .);j[) 



(3.85) 
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where C(s, r) = 32^ 



VT- 



The idea for obtaining exponential decay is to use the terms E(| (Cl"^ — 
^)ksiY) *o smallness and the terms E(|(C|"^ — to repeat the 

process. Thus, using the Lemma ESI we get that for any < 1, there 
exists a fixed constant m > such that, for any si, \si — {k + m)\ < 2, 
we have 

4^ ■ 64 • r) ■ E (I (C;-) - ^)-\ [) < /? (3.86) 

We can now repeat the same procedure for each term E(|(Cj"'' — 
z)~^iY) and we gain one more coefficient (3. At each step, we move 
(m + 3) spots to the right from k to /. We can repeat this procedure 
[m^] times and we get 

E (^|(C(") - z)fc/f) < (3.87) 

which immediately gives ()2.6|1 . □ 



4. The Localized Structure of the Eigenfunctions. 

In this section we will study the eigenfunctions of the random CMV 
matrices considered in (j2.ip . We will prove that, with probability 1, 
each eigenfunction of these matrices will be exponentially localized 
about a certain point, called the center of localization. We will fol- 
low ideas from del Rio et al. [5]. 

Theorem 12.11 will give that, for any z G c?D, any integer n and any 

se(0,i), 

E(|FH(^,CW)f) <Ce-^l^-^l (4.1) 

Aizenman's theorem for CMV matrices (see Simon ^j) shows that 
f|4.ip implies that for some positive constants Co and Dq depending on 
s, we have 

E fsup I (4, (Ci'^V'^z) \] < Co e~^"l'-^l (4.2) 

This will allow us to conclude that the eigenfunctions of the CMV 
matrix are exponentially localized. The first step will be: 

Lemma 4.1. For almost every a E Q, there exists a constant Da > 
such that for any n, any k, I, with 1 < k,l < n, we have 

snp\{6k,{Cj^'>y6i)\<Da{l + n)'e-'''>\'-'\ (4.3) 
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Proof. From ()4.2j) we get that 

/ f sup I {6k, (C("))^5^) dPia) < Co e'^^l'^-^l (4.4) 
and therefore there exists a constant Ci > such that 

(4,5) 

It is clear that for any k, I, with 1 < k,l < n, the function 
is integrable. 

Hence, for almost every a & fl, there exists a constant > such 
that for any n, k, I, with 1 < k,l < n, 

sup|(5,,(C("))^5z)| < Z^Jl + nf e-^^l'^-'l (4.7) 

□ 

A useful version of the previous lemma is: 

Lemma 4.2. For almost every a & Q, there exists a constant > 

such that for any n, any k, I, with 1 < k,l < n, and \k — l\ > ln(n + 
1), we have 

sup\{S,,{Cj^^ySi)\<C^e~^^'-'^ (4.8) 



Proof. It is clear that for any n, k, I, with 1 < k,l < n and \k — l\ > 
f ln(n + l), 

1 ,^\k-l\ 



(l + n2)(l + A;2)(l + /2) 



e 2 I'^-'l > 1 (4.9) 



In particular, for any n, k, I with \k — l\ > \n{n + 1), the function 



n3a ^ fsup|(5fc,(C('^))^5i)|)e^l'=-'l (4.10) 
\j& / 

is integrable, so it is finite for almost every a. 

Hence for almost every a & Q, there exists a constant Ca > such 

that for any k,l, \k — l\ > j^^ ln(n + 1), 

sup|(5,,(Ci"))^50| < C^e-^l'^-'l (4.11) 

□ 



24 MIHAI STOICIU 

Proof of Theorem \2.^ Let us start with a CMV matrix C*-"-* = Ci"'' 
corresponding to the Verblunsky coefficients ao, ai, . . . , an-2, As 
mentioned before, the spectrum of C^"^ is simple. Let e*^° be an eigen- 
value of the matrix Ci"'' and Lp^a'^ a corresponding eigenfunction. 
We see that, on the unit circle, the sequence of functions 

M 

^ ^ 2M + 1 ^ ^ ' 

j=-M 

is uniformly bounded (by 1) and converge pointwise (as M —>■ oo) to 
the characteristic function of the point e*^". Let P^e^Oa} = X{eiOa}{C^^)- 



By Lemma [4.21 we have, for any k, I, with \k — l\ > ln{n + 1) 

1 



Do 

M 

|(5.,/m(c("V0 



2M 



(4.13) 



i=-M 
M 

" |fc-Z| 



j=-M 



(4.14) 



where for the last inequality we used ()4.1|) . 

By taking M — > oo in the previous inequality, we get 



|(4,P{e»«.}MI < C^e"^\^-'\ (4.15) 
and therefore 

ki"Hfc)^i"HO| < C^e~^\^-'\ (4.16) 
We can now pick as the center of localization the smallest integer 
m{'~pi^^) such that 

I^L"H^(^L"^))l=max|^(r)(m)| (4.17) 

m 

We clearly have \p^a\m{Lp^a'^))\ > -^t^^- 



Using the inequality ()4.1f)j) with k = m and / = m{ipa ) we get, for 
any m with |m — ■m{ip^a^)\ > \n{n + 1), 

\'p'a\m)\< C7^e-^l"(^-"')l V^TTT (4.18) 

Since for large n, e~^''^~'' y/n + 1 < e~^l^~'l for any k,l, \k — l\ > 
ln(n + 1), we get the desired conclusion (we can take D2 = j^). □ 
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For any eigenfunction , the point m{(pa ) is called its center of 
localization. The eigenfunction is concentrated (has its large values) 
near the point m{(f^a^) and is tiny at sites that are far from m{(f^a^). 
This structure of the eigenfunctions will allow us to prove a decoupling 
property of the CMV matrix. 

Note that we used Lemma 14.21 in the proof of Theorem 12.21 We can 
get a stronger result by using Lemma HH] (we replace ()4.11|) by ()4.7|) ). 
Thus, for any n and any m < n, we have 

\v'a\^)\< /^a(l + n)''e-^l"(^""^)l V^TTT (4.19) 
where m(<^i"'') is the center of localization of the eigenfunction ip^a^ . 



5. Decoupling the Point Process 

We will now show that the distribution of the eigenvalues of the CMV 
matrix C^"-* can be approximated (as n — > oo) by the distribution of 
the eigenvalues of another matrix CMV matrix C^"-*, which decouples 
into the direct sum of smaller matrices. 

As explained in Section^ for the CMV matrix C*^"-* obtained with the 
Verblunsky coefficients a = (ao, ai, . . . , an-i) G we consider C^^^ the 
CMV matrix obtained from the same Verblunsky coefficients with the 
additional restrictions ar^i = e*''^ Q^2f^l ~ • • • ' ^n-i = e*'''['""i, 

[ In ri J [ In n J 

where e*''^, e*''^, . . . , e*'^['""i are independent random points uniformly 
distributed on the unit circle. The matrix C'-"^ decouples into the di- 
rect sum of approximately [Inn] unitary matrices C["'\ C'^\ . . . , Cp^^n]- 
Since we are interested in the asymptotic distribution of the eigenval- 
ues, it will be enough to study the distribution (as n — > oo) of the 
eigenvalues of the matrices C*-^-* of size = [Inn] [j^]. Note that in 
this situation the corresponding truncated matrix C'-^-' will decouple 
into the direct sum of exactly [Inn] identical blocks of size [j^]- 
We will begin by comparing the matrices C*-^^ and C^^\ 

Lemma 5.1. For N = [Inn] [j^], the matrix C^^^ —C^^^ has at most 
4[lnn] nonzero rows. 

Proof. In our analysis, we will start counting the rows of the CMV 
matrix with row 0. A simple inspection of the CMV matrix shows 
that for even Verblunsky coefficients a2k, only the rows 2k and 2k + 1 
depend on a2k- For odd Verblunsky coefficients a2k+i, only the rows 
2k, 2k + 1,2k + 2, 2k + 3 depend on a2k+i- 
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Since in order to obtain the matrix C^^^ from C*-^-* we modify [Inn] 
Verblunsky coefficients ar^i, a2r^], . . . , «nn„ir_2_i, we immediately 

see that at most 4[ln?T,] rows of C^^^ are modified. 

Therefore C-^-* — C*-^-* has at most 4[lnn] nonzero rows (and, by the 
same argument, at most 4 columns around each point where the matrix 
C^^^ decouples). □ 

Since we are interested in the points situated near the places where 
the matrix C^^^ decouples, a useful notation will be 

Sn{K) = S^^\K) U S^^\K) U • • ■ U (5.1) 

where S'-''\K) is a set of K integers centered at k [j^] (e.g., for 
K = 2p, S(^\K) = {k[^]-p + l,k[^]-p + 2,... T[^] +p}). 
Using this notation, we also have 

Sn{1) = \ r— ,2 — ,...,[lnn] — 5.2 

Consider the intervals lN,k, 1 < k < m, of size near the point 

on the unit circle (for example /tv^^ = (e**-""'""iv e*'-""'"^-')), where 
cti < bi < a2 < b2 < ■ ■ ■ < a^n < bm- We will denote by A/'iv(/) the 
number of eigenvalues of C^'^^ situated in the interval /, and by J\fN{I) 
the number of eigenvalues of C-^^ situated in /. We will prove that, for 
large A^, J^N{lN,k) can be approximated by AfNilN,k), that is, for any 
integers ki, k2, ■ ■ ■ , km > 0, we have, for N oo, 

I P(Arjv(/jv,l) = fcl, Ar^(/iV,2) = k2,.. ■,AfNilN,m) = k^) 
- ¥{Mn{In,i) = ki, Mn{In,2) = k2,.. .,UN{lN.m) = /^m) | ^ 

(5.3) 

Since, by the results in Section 01 the eigenfunctions of the matrix 
are exponentially localized (supported on a set of size 2T[ln(n + 
1)], where, from now on, T = -^), some of them will have the cen- 
ter of localization near Sn{^) (the set of points where the matrix 
C^^^ decouples) and others will have centers of localization away from 
this set (i.e., because of exponential localization, inside an interval 
{k [^].{k + l) [^])). 

Roughly speaking, each eigenfunction of the second type will produce 
an "almost" eigenfunction for one of the blocks of the decoupled matrix 
C^^\ These eigenfunctions will allow us to compare N'N{lN,k) and 

A/7v(/Ar,A:)- 

We see that any eigenfunction with the center of localization outside 
the set S'Ar(4T[lnn]) will be tiny on the set S]^{1). Therefore, if we want 
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to estimate the number of eigenfunctions that are supported close to 
Sn{1), it will be enough to analyze the number bN^a, where bj^^a = 
number of eigenfunctions of with the center of localization inside 
5'7v(4T[lnn]) (we will call these eigenfunctions "bad eigenfunctions"). 
We will now prove that the number bj^ ^ is small compared to A^. 

A technical complication is generated by the fact that in the expo- 
nential localization of eigenfunctions given by (j4.H|l . the constant Da 
depends on a E Q. We define 

MK=Uen, sup|(4,(C(^)WI <^(l + iVfe-^°l^-'lj (5.4) 
I j& ) 

Note that for any > 0, the set TW^^ C is invariant under ro- 
tation. Also, we can immediately see that the sets Aix grow with K 
and 

lim F{Mk) = 1 (5.5) 

We will be able to control the number of "bad eigenfunctions" for 
a & Aix using the following lemma: 

Lemma 5.2. For any K > and any a E M.k, there exists a constant 
Ck > such that 

bN,a<CK{Hl + N)f (5.6) 

Proof. For any K > 0, any a E Aix, and any eigenfunction which 
is exponentially localized about a point m{ip^), we have, using ()4.19|) . 

\ipl^\m)\< ire-^l"(^"''')l (l + Ar)*^v^rTiV (5.7) 

Therefore for any m such that \m-m{ip^)\ > \^\n{l + N) 
have 



\m-m{vS)\>[^ln{l+N)\ 

<(l + Ar)-ir^-^ (5.8) 

Therefore, for any fixed K and s, we can find an A^o = No{k, s) such 
that for any N > Nq, 

E k'I'\rn)\'>l (5.9) 

|m-m«)|<[^ln(l+iV)] 

We will consider eigenfunctions (f^ with the center of localization in 
S'at (4T[ln A^]). For a fixed a G M.k, we denote the number of these 
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eigenfunctions by b^^a- We denote by {4'i,4'2, ■ ■ ■ yi^bM A the set of 
these eigenfunctions. Since the spectrum of C-^^ is simple, this is an 
orthonormal set. 

Therefore if we denote by card(y4) the number of elements of the set 
A, we get 

< cardj^ (^4r[lnAr] + ^\n{l + N) 
14\ 2 



<[4^T+—j (ln(l + iV)) 
Also, from dSH), for any N > No{K, s), 

E E (5.10) 

mes(4T[lnAr] + [-i^ln(l+Af)]) 

Therefore, for any i^' > and any a G M.Ki we have, for N > 
No{K,s), 

bN,a < 2 (^4T + -^^ (ln(l + N)Y (5.11) 

and we can now conclude ()5.6|) . □ 

Lemma shows that for any K > 0, the number of "bad eigenfunc- 
tions" corresponding to a G TWii- is of the order (InA^)^ (hence small 
compared to N). 

Since the distributions for our Verblunsky coefficients are taken to 
be rotationally invariant, the distribution of the eigenvalues is rota- 
tionally invariant. Therefore, for any interval In of size on the unit 
circle, and for any fixed set M-k C Q, the expected number of "bad 
eigenfunctions" corresponding to eigenvalues in 1^ is of size ^^^^^ ■ We 
then get that the probability of the event "there are bad eigenfunctions 
corresponding to eigenvalues in the interval In" converges to 0. This 
fact will allow us to prove 

Lemma 5.3. For any K > 0, any disjoint intervals In,i, In,2i ■ ■ ■ ■, lN,m 
{each one of size jj and situated near the point e*") and any positive 
integers ki, k2, ■ ■ ■ , km, we have 

I F({Ar^(J^,i) = ki, AfNilN,2) =k2,.. . , A/'7v(/7V,m) = km} H Mk) 
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- F{{J^n{In,i) = h, AGv(/iV,2) =k2, ... AfNilN^m) = km} H Mk) \ ^ 

(5.12) 

as ^ oo. 

Proof. We will work with a & A4k- We first observe that any "good 
eigenfunction" (i.e., an eigenfunction with the center of localization 
outside Sn (4T[lnA^])) is tiny on 5'Ar(l). 

Indeed, from ()4.19j) . for any eigenfunction y^i^'' with the center of 
localization m{ip"a^) and for any m with \'m—m{if"a'')\ > -^[ln(A^+l)], 

\ip\^\m)\ < i^e-^l"(^"'")l(l + Nf^/lTN (5.13) 

In particular, if the center of localization of tp^a'^ is outside (4T[ln A^]), 
then for all m G S^iX), we have 

\^\^\m)\ <K(l + iV)-2 (5.14) 

We will use the fact that if is a normal matrix, 2:0 G C, e > 0, and 
ip is a. unit vector with 

II (AT - zo)v? II <e (5.15) 

then A^ has an eigenvalue in {2; | |z — zqI < 

For any "good eigenfunction" 'p>^a\ we have C^^(p^^ = and there- 
fore, using Lemma f5.H 

IICfVL'^^ll < 2K[ln Ar](l + N)-^ (5.16) 
Therefore, for any interval In of size j^, we have 

^fNiIN)<^fNi^N) (5.17) 

where 1^ is the interval 1^ augmented by 2i^[ln A^](l + A^)~^. 
Since 2_ft'[ln A^](l + A^)~^ = o(-^), we can now conclude that 

P((A/iv(/iv) < J^NilN)) n Mk) -^1 as n^oo (5.18) 

We can use the same argument (starting from the eigenfunctions of 
which are also exponentially localized) to show that 

P((7Vjv(/jv) > J^NilN)) n Mk) ^1 as n ^ oo (5.19) 

so we can now conclude that 

P((Ar^(/iv) = Ar^(/iv)) n Mk) ^1 as n ^ oo (5.20) 

Instead of one interval In, we can take m intervals In,i, In, 2 • • • , lN,m 
so we get (jSUll). □ 
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Proof of Theorem, \2.'A Lemma 15.31 shows that for any > 0, the dis- 
tribution of the eigenvalues of the matrix C^^^ can be approximated by 
the distribution of the eigenvalues of the matrix C'-^^ when we restrict 
to the set J^k C VL. Since by fIS.Sp the sets M.k grow with K and 
limi5:^oo IP(A^k) = 1, we get the desired result. □ 



6. Estimating the Probability of Having Two or More 
Eigenvalues in an Interval 

The results from the previous section show that the local distribu- 
tion of the eigenvalues of the matrix C*^^-* can be approximated by 
the direct sum of the local distribution of [In n] matrices of size [j^] , 

Cj^'' , , • • • , C,\^n] ■ These matrices are decoupled and depend on in- 
dependent sets of Verblunsky coefficients; hence they are independent. 

For a fixed point e'^° e (9©, and an interval 1^ = (e*(^o+^\ e*(^o+"^)), 
we will now want to control the probability of the event "C'-^^ has k 
eigenvalues in 1^." We will analyze the distribution of the eigenvalues 
of the direct sum of the matrices Cj^-* , €2^^ , • • • , C^^^j . We will prove 

that, as n 00, each of the decoupled matrices C^^'' can contribute 
(up to a negligible error) with at most one eigenvalue in the interval 
In- 

For any nonnegative integer m, denote by A{m,C,I) the event 

A{m, C, I) = "C has at least m eigenvalues in the interval /" (6.1) 

and by i?(m, C, /) the event 

B{m, C, I) = "C has exactly m eigenvalues in the interval /" (6.2) 

In order to simplify future notations, for any point e*^ G 9D, we also 
define the event M(e*^) to be 

M{e'^) = "e*^ is an eigenvalue of C^^^" (6.3) 

We can begin by observing that the eigenvalues of the matrix C^^-* 
are the zeros of the A^-th paraorthogonal polynomial (see (jl.2|) ): 

^N{z,diJ,f3) = z^N~i{z,dfi) - ^^%^^{z,di2) (6.4) 

Therefore we can consider the complex function 

- (6.5) 

which has the property that $Ar(e*^) = if and only if i?Ar(e*^) = 1. 
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By writing the polynomials ^n-i and products of their 

zeros, we can see that the function Bj^ is a Blaschke product. 
Let r^jv : [0, 2tt) ^ M be a continuous function such that 

B^ie'^) = e^''^^'^^ (6.6) 

(we will only be interested in the values of the function rj^ near a fixed 
point e*^° G 83). Note that for any fixed 9 G 83, we have that ri{6) is a 
random variable depending on a = (ao, C(i, ■ ■ ■ , o:n-2, o:n-i = P) E fl. 

We will now study the properties of the random variable rii\i{6) = 
riN{e, ao, aN-2, P)- Thus 

Lemma 6.1. For any 6i and 62, the random variables ^^{61) and 
T]N{d2) are independent. Also for any fixed value w e 

^(^((^i) VN{e2) = w] =N (6.7) 



Proof. The equation ()6.5|1 gives 

VN{0)=^ + T{e) (6.8) 

where e''"' = (3 and e*"^*^^^ = ^ ^?^~(e»d) • Since the distribution of each of 

the random variables ao, oti-, ■ ■ ■ , otN-2 and (3 is rotationally invariant, 
for any 9 G [0,27r), 7 and t{9) are random variables uniformly dis- 
tributed. Also, it is immediate that 7 and t{9) are independent. Since 
7 does not depend on 9, for any fixed ^1, 6*2 G [0, 27r), we have that the 
random variables ^-(6'i) and r]N{92) are independent. 

We see now that for any Blaschke factor Ba{z) = fz^, we can define 
a real-valued function rja on (9D such that 

A straightforward computation gives 

Since B^q is a Blaschke product, we now get that for any fixed a G fi, 
has a constant sign (positive). This implies that the function rjN 
is strictly increasing. The function Bn{z) is analytic and has exactly 
N zeros in D and therefore we get, using the argument principle, that 

J^'^^{9)d9 = 27rN (6.11) 

Note that does not depend on /? (it depends only on oq, cti, . . . , aAr_ 
Also, using the same argument as in Lemma f3.H we have that for any 
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angles 9 and 

(6.12) 

where fj is the function rj that corresponds to the Verblunsky coefficients 

afc,^ = e-^(^+i)^afc = 0, 1, . . . , (iV - 2) (6.13) 

Since the distribution of ao, ai, . • • , aiv-2 is rotationally invariant, 
we get from (jO^ that the function ^ ^ E {^{0)) is constant. 

Taking expectations and using Fubini's theorem (as we also did in 
Lemma EIH), we get, for any angle 6'o, 

(6.14) 

and therefore 

E (^(^^o)) = N (6.15) 

Since for any 6^1,^2 ^ [0,27r), we have that ^-(^1) and ?7Ar(^2) are 
independent, ()6.15|) implies that for any fixed value w G M, 

Ef^(^^i) VN{e2) = w]=N (6.16) 
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□ 



We will now control the probability of having at least two eigenvalues 
in /jv conditioned by the event that we already have an eigenvalue at 
one fixed point e*^^ G In- This will be shown in the following lemma: 

Lemma 6.2. WithC^^\lM, and the events A{m,C, I) and M{e'^) de- 
fined before, and for any e*^^ G In, we have 

P(A(2,C(^),Jjv) I M{e'^)) <{b-a) (6.17) 

Proof. Using the fact that the function 9 ^ K {^§^{9)) is constant and 
the relation ()6.16|) . we get that 

= w\ =2n{b-a) (6.18) 




^N{en = ^ BN{e'') = 1 ^ t]n{9) = (mod 27r) (6.19) 

Therefore if the event A{2,C^^\ In) takes place (i.e., if the polyno- 
mial $Ar vanishes at least twice in the interval In), then the function 
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?7Ar changes by at least 2ti in the interval In, and therefore we have 
that 



, r,Ad)de >27r (6.20) 

whenever the event A{2,C^'^\ In) takes place. 

For any 6i G /at we have, using the independence of the random 
variables ^^{Oi) and r]N{(^2) for the first inequality and Chebyshev's 
inequality for the second inequality, 



P (a{2,C^^\In) M{e'^')) <f( [ ^{0) dO > 27r M{e'^'' 

2tt \Jeo+^ 96 

(6.21) 

The previous formula shows that we can control the probability of 
having more than two eigenvalues in the interval In conditioned by the 
event that a fixed e*^^ is an eigenvalue. We now obtain, using ()6.18p 
with w = 27rm, m G Z, 

¥{A{2,C^''\In) I M{e''')) < {b - a) (6.22) 

□ 

We can now control the probability of having two or more eigenvalues 
in In- 

Theorem 6.3. With C^^\ In, and the event A{m, C, /) defined before, 
we have 

P(A(2,C(^),J^)) <^^^ (6.23) 
Proof. For any positive integer k, we have 

P(fi(A;,C(^),/^)) = - / P(i?(A;,C(^),/^) | M(e^^)) N duNiO) 

(6.24) 

(where the measure un is the density of eigenvalues). 

Note that the factor ^ appears because the selected point e*^ where 
we take the conditional probability can be any one of the k points. 

We will now use the fact that the distribution of the Verblunsky 
coefficients is rotationally invariant and therefore for any we have 



34 MIHAI STOICIU 

dfN = where ^ is the normahzed Lebesgue measure on the unit 
circle. 

Since for any > 2 we have | < |, we get that for any integer k > 2 
and for large A^, 

N _ , ... .... dd 



P(i?(A:,C(^),/^)) < - / P(i?(fc,C(^),/;v) I M(e^^)) 



2n 

(6.25) 



and therefore, 

^ _|_ 27vb 

P(A(2,C(^),/^)) < ^ / ° P(A(2,C(^),/^) I M(e^^)) f 

2 Jeo+2fi ' ^ 27r 

(6.26) 

Using Lemma (6. 2^ we get 

P(^(2.C«M)<^ii9^((,-a)=(^^ (6,27) 

□ 

Theorem 6.4. Mt/i C(^), C(^\ C^^\ . . . , Cj^^p J^v, and the event A{m, C, /) 
defined before, we have, for any k, 1 < k < [Inn], 

P(^A(2,Cf\/jv)) =0(([lnn])~') as n ^ oo (6.28) 

Proof. We will use the previous theorems for the CMV matrix Cl^\ 
Recall that = [Inn] [j^]- Since this matrix has [j^] eigenvalues, 
we can use the proof of Lemma f6. 21 to obtain that for any e*^ G In, 

1 27r(6-a) 



p(.(2,Cr,«|M,e.«,).^^[^]^^ (-) 
The proof of Theorem 16.31 now gives 

and hence ()6.28p follows. □ 

This theorem shows that as — *• cxo, any of the decoupled matrices 
contributes with at most one eigenvalue in each interval of size 

7. Proof of the Main Theorem 

We will now use the results of SectionsEHHEl andElto conclude that 
the statistical distribution of the zeros of the random paraorthogonal 
polynomials is Poisson. 
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Proof of Theorem M.ll It is enough to study the statistical distribution 
of the zeros of polynomials of degree N = [In n] [ • These zeros 
are exactly the eigenvalues of the CMV matrix C^^\ so, by the results 
in Section El the distribution of these zeros can be approximated by 
the distribution of the direct sum of the eigenvalues of [Inn] matrices 

MN) MN) MN) 
l-l 51-2 ' • • • ' '-[Inn]- 

In Section lUl (Theorem I6.4j) . we showed that the probability that any 
of the matrices C[^^ , Cg^"* , . . . , C^l^^^^ contributes with two or more eigen- 
values in each interval of size situated near a fixed point e*^ G dlD) is 
of order 0([ln Since the matrices C['^\C2'^\ . . . , C^^.^^ are identi- 

cally distributed and independent, we immediately get that the proba- 
bility that the direct sum of these matrices has two or more eigenvalues 
in an interval of size situated near e*^ is [Inn] 0([lnn]~^) and there- 
fore converges to as — * 00. 

We can now conclude that as n ^ 00, the local distribution of the 
eigenvalues converges to a Poisson process with intensity measure n ^ 
using a standard technique in probability theory. We first fix an interval 
In = {e''^^°~^^\ e*(^''+"^)) near the point e^^° (as before, we take = 
[Inn] [i^]). Let us consider [Inn] random variables Xi,X2, . . . ,X[in„] 
where = number of the eigenvalues of the matrix C^.^'* situated in 
the interval In and let Sn{lN) = Xi + X2 + ■ ■ ■ + X[\^n]- Note that 
Sn{lN) = the number of eigenvalues of the matrix C^^^ situated in the 
interval In- We want to prove that 

hm P(^.(/^) = k) = e-(^-»)i^_^ (7.1) 

n— >oo k\ 

Theorem 16.41 shows that we can assume without loss of generality 
that for any k,l < k < [Inn], we have X^ G {0, 1}. Also, because of 
rotation invariance, we can assume, for large n, 

= 1) = ^ (7.2) 
[m n] 

P(X, = 0) = 1 - (7.3) 
[mn] 

The random variable Sn{lN) can now be viewed as the sum of [Inn] 
Bernoulli trials, each with the probability of success and 
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which converges to e where A = [Innj^p^ = (6 — a). Therefore 
we get (jZH). 

Since for any disjoint intervals lN,k, 1 < ^ < [Inn] situated near e*^°, 
the random variables Sn{lN,k) are independent, ()7.H) will now give p.6p 
and therefore the proof of the main theorem is complete. □ 



8. Remarks 

1. We should emphasize the fact that the distribution of our random 
Verblunsky coefficients is rotationally invariant. This assumption is 
used in several places and seems vital for our approach. It is not 
clear how (or whether) the approach presented here can be extended 
to distributions that are not rotationally invariant. 

2. In this paper, we study the statistical distribution of the zeros 
of paraorthogonal polynomials. It would be interesting to understand 
the statistical distribution of the zeros of orthogonal polynomials. A 
generic plot of the zeros of paraorthogonal polynomials versus the zeros 
of orthogonal polynomials is 
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In this Mathematica plot, the points represent the zeros of paraorthog- 
onal polynomials obtained by randomly choosing ao, ai, . . . , agg from 
the uniform distribution on D{0, |) and ajQ from the uniform distri- 
bution on (9D. The crosses represent the zeros of the orthogonal poly- 
nomials obtained from the same Uq, ai, . . . , agg and an ajo randomly 
chosen from the uniform distribution on D{0,^). 
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We observe that, with the exception of a few points (corresponding 
probably to "bad eigenfunctions" ) , the zeros of paraorthogonal polyno- 
mials and those of orthogonal polynomials are very close. We conjec- 
ture that these zeros are pairwise exponentially close with the exception 
of 0((lniV)^) of them. We expect that the distribution of the argu- 
ments of the zeros of orthogonal polynomials on the unit circle is also 
Poisson. 

3. We would also like to mention the related work of Bourget, How- 
land, and Joye [3^ and of Joye [1^], [IE]- In these papers, the authors 
analyze the spectral properties of a class of five-diagonal random uni- 
tary matrices similar to the CMV matrices (with the difference that 
it contains an extra random parameter). In JE] (a preprint which 
appeared as this work was being completed), the author considers a 
subclass of the class of |5j that does not overlap with the orthogonal 
polynomials on the unit circle and proves Aizenman-Molchanov bounds 
similar to the ones we have in Section El 

4. The results presented in our paper were announced by Simon 
in |23j, where he describes the distribution of the zeros of orthogonal 
polynomials on the unit circle in two distinct (and, in a certain way, op- 
posite) situations. In the first case, of random Verblunsky coefficients, 
our paper shows that there is no local correlation between the zeros 
(Poisson behavior). The second case consists of Verblusky coefficients 
given by the formula a„ = Cfe" + ©((feA)") (i.e., converges to 
a constant C sufficiently fast). In this case it is shown in that 
the zeros of the orthogonal polynomials are equally spaced on the cir- 
cle of radius 6, which is equivalent to saying that the angular distance 
between nearby zeros is 2TT/n ("clock" behavior). 
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